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In this paper, we investigate the non-Markovian dynamics of a model to mimic the excitation en- 
ergy transfer (EET) between chromophores in photosynthesis systems. The numerical path integral 
method is used. This method includes the non-Markovian effects of the environmental affects and 
it does not need the perturbation approximation in solving the dynamics of systems of interest. It 
implies that the coherence helps the EET between chromophores through lasting the transfer time 
rather than enhances the transfer rate of the EET. In particular, the non-Markovian environment 
greatly increase the efficiency of the EET in the photosynthesis systems. 
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I. INTRODUCTION 

Photosynthesis provides almost all the energy for life 
on Earth. Therefore, many interests have been attracted 
on this topic in past decades not only in biology but 
also in chemistry and physics [l[ . The typical antenna of 
harvesting photons in photosynthesis systems are protein 
complexes that hold pigments [chlorophyll (Chi), bacte- 
riochtorophyll (BChl), phycobilin, carotenoid, bacterio- 
pheophytin (BPhy), etc.]. For example, in a green sulfur 
bacteria, the Fenna-Matthews-Olson (FMO) protein is a 
trimer made of identical subunits, each of which contains 
seven BChl molecules @. Another purple bacterium, Rh. 
sphaeroides @ includes a special pair Chls (P) in the cen- 
ter, accessory Bchls (B^ and Bju) on the each side of P, 
and BPhys (Rl and Hm) next to the each BChls. 

The photosynthesis starts with the absorption of pho- 
ton of sunlight by the light-harvesting pigments, followed 
by the excitation energy transfer (EET) to the reaction 
center, where charge separation is initiated and physi- 
cal energy is converted into chemical energy Q. At low 
light intensities, the quantum efficiency of the transfer 
of energy is very high, with a near unity quantum yield. 
Although this process has been investigated extensively 
the details and the mechanism of the energy transfer has 
not been fully understood today. 

There is a hypothesis that the high efficiency of near 
unity EET attributes to coherent superpositions between 
excitation states that belong to different chromophores 
of the photosynthesis systems. It is deduced from the 
following analogous facts. Engel et al. 5j] found that in 
FMO protein the coherence between the two excitation 
states clearly lasts for a longer time which is similar to 
the time scale of the EET. Lee et al.[6| also found that in 
Rh. sphaeroides the coherence of the excitations created 
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from laser pulses on adjacent Br, and H^ can prolong for 
longer time which is long enough for the coherent trans- 
fer of energies between chromophores B l and H l ■ These 
kinds of experimental results imply that electronic excita- 
tions may move coherently through these photosynthesis 
systems rather than by incoherent hopping motion as has 
usually been assumed. 

In the aspect of theory, more than fifty years earlier, 
ForsterQ suggested a formula of the efficiency of the 
resonance EET. This theory give a picture that excita- 
tion energies are transferred in the photosynthesis sys- 
tems through incoherent hopping from the higher levels 
to the lower levels. However, this theory only describes 
the cases that the couplings between chromophores are 
weak. The master equation of Redfield form Q has also 
been used in the investigations of the EET. However, 
this method of the original form is limited because it 
derived from the Markovian and perturbation approxi- 
mations. A modified Redfield theory, first suggested by 
Zhang and Fleming Q, and later elaborated by Yang 
and Fleming [Io| . has then been used in the investiga- 
tions of the EET. It partly overcomes the shortcomings 
of the conventional one. But, it cannot take into ac- 
count any processes related to quantum coherence and 
decoherence. Jang et al. [ll[ developed the F orster- 
Dexter theory for the EET. In their method the weak 
couplings between chromophores are not needed, but it 
asks the couplings between the system and the environ- 
ment be weak. Ishizaki and Fleming [l2j recently used 
the reduced hierarchy equation approach dealing with the 
quantum coherent and incoherent dynamics of the EET. 
It avoids the usage of perturbative truncation, and can 
describe the quantum coherent wavelike motion and in- 
coherent hopping in the same framework. Mohseni et 
al. |13l | developed a quantum walk approach, based on a 
quantum trajectory picture in Born-Markovian and sec- 
ular approximations, as a natural framework for incorpo- 
rating quantum dynamical effects in the EET. It is shown 
that the interplay between the coherent dynamics of the 
system and the incoherent action of the environment can 
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lead to significantly greater transport efficiency than the 
coherent dynamics on its own. Similar results have also 
been obtained independently by Plenio et al.[Ji|. Much 
other effort bas been contributed to the fields in last years 

EMI. 

However, in the existing investigations as being listed 
above, two important aspects paid not enough attention 
to. One is how the non-Markovian [20|, |2l| effects influ- 
ences the dynamics of the EET, and the other is how to 
describe the evolution rate of the EET. In this paper, we 
shall use the numerical path integral method developed 
by Makri's group investigating the EET, where we will 
pay particular attention to above two problems. 

This article is organized as follows. In Sec. II, we 
introduce the model and compare the numerical path in- 
tegral method to the master equation method through 
calculating the time evolution of the population of sites. 
In Sec. Ill we investigate the evolutions of the coherent 
term of the reduced density matrix and the transfer rate 
of the EET in many kinds of conditions. It will be shown 
that the coherence helps the energy transfer in the model 
system through lasting the transfer time rather than in- 
creasing the transfer rate of the EET, and as we consider 
the non-Markovian effects the energy transfer will be in- 
creased greatly. Through the model study we can obtain 
some important implications for energy transfer in pho- 
tosynthesis systems. Finally, Sec. IV is devoted to the 
concluding remarks. 



II. MODEL AND TIME EVOLUTION OF THE 
POPULATION OF SITES 

The Hamiltonian of the EET model can be represented 
by Frenkel exciton Hamiltonian [l(| E3 • 



H fnf = H el + H ph + H re + H el - ph , 



(1) 



where 



Hel = Eb') £ °01+ A (l 1 )( 2 l + I 2 )( 1 l)' ( 2 ) 
i=i 

2 
2 

H re = £|j)A,(j|, (4) 
i=i 



i=i j'=i 



H el-ph = J2 Hf- ph = Yi 



(5) 



Here, |j) represents the state where only the jth site is 
in its excitation electronic state \tfij e ) and another is in 
its ground electronic state \fkg), i-e., |j) = \fj e ) Wkg)- 
While we define the ground state |0) = \<£i g ) \y>2g)- 
Eq. ([2]) is the electronic Hamiltonian in which the Hamil- 
tonian of the ground electronic state is set to be zero. 



is the excitation electronic energy of the jth site in 
the absence of phonon. A is the electronic coupling be- 
tween the two sites, which responsible for the EET be- 
tween the individual sites. H ph = £^ hu^ (p 2 + q"^j /2 

is the phonon Hamiltonian associated with the jth site, 
where q^, and u>£ are the dimensionless coordinate, 
conjugate momentum, and frequency of the £th phonon 
mode, respectively. Xj = ~^2^fiio^d 2 ^/2 is the reorganiza- 
tion energy of the jth site, where dj£ is the dimensionless 
displacement of the equilibrium configuration of the £th 
phonon mode between the ground and excitation states 
of the jth site. H el ~ ph is the coupling Hamiltonian be- 
tween the jth site and phonon modes and the Vj and 
Uj are defined as Vj — |j) (j| and Uj — hoj^dj^. 
Here, we assume the phonon modes associated with one 
site are uncorrelated with those of another site. If wc 
set Ai = A2 = A, namely, suppose that the reorganiza- 
tion energy of the 1-th and 2-th sites are equal, and reset 
the zero point energy at (e" — s\ + u 2 — «i)/2, and set 
C£ = (^2{ — d\£) /2, we can obtain the total Hamiltonian 
presented with Pauli matrix as 



H 



( -<j z + Aa x - a z c ( Uuj ( q^ 



(6) 



where e = £2~ e i- Eq-® is the standard form of the spin- 
boson model. Our following investigations start from 
this Hamiltonian. Thus, we can investigate the model 
with numerical path integral scheme which is developed 
by Makri group (22j |. and was widely used by Thorwart 
group [23| and many other groups [24| in last years. 

In order to compare numerical path integral scheme 
to master equations approach of Redfield form we firstly 
investigate the time evolution of the population of sites 
for the model. The spectral distribution function is also 
employed as the Drude-Lorentz density (the over-damped 
Brownian oscillator model): 



J{uj) = 2A- 



LU", 



(7) 



as Ref.[l2j. The parameters are taken as the same val- 
ues as Ref.[l2;|, namely, e = 100 cm -1 , A = 100 cm -1 , 
7 = 53.08 cm" 1 (7- 1 = 100 fs), and T = 300 k. To 
implement the numerical path integral scheme, we use 
the iterative tensor multiplication (ITM) algorithm de- 
rived from the the quasiadiabatic propagator path inte- 
gral (QUAPI) technique. We use the time step St = 10 fs 
which is shorter than the correlation time of the bath and 
the characteristic time of the two-level system. In order 
to include as much non-Markovian effect of the bath as 
possible in the ITM algorithm, one should choose Sk max 
as large as possible so that the dominant non-Markovian 
effect of the bath is included. Here, Sk max is the maxi- 
mal number of time steps, and the 5£; max <5i is the included 
memory time in the calculations. In another word, the 
memory effects of the bath within time Sk max St are con- 
sidered in the algorithm. If set <5fc ma x = 0, this method 
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FIG. 1: (Color online) Time evolutions of the population of 
site 1, obtained from master equations of Redfield form (black 
line connects black squares), and from ITM algorithm based 
on QUAPI of 5kmax = (blue dashed line connects blue 
stars), and 8k max = 3 (red dotted line connects red circles) 
for different A. Here, e = 100 cm -1 , A = 100 cm" 1 , 7 = 53.08 
cm -1 (7 _1 = 100 fs), and T — 300 k. Initially, the system is 
populated on the site 1 with state |-0o) = !)• 



then reduce to the Markovian approximation. But, the 
Markovian results are not same to the ones from master 
equations of Redfield form, because the latter is not only 
Markovian but also perturbative. The time evolutions of 
the population of site 1 (suppose at the initial time the 
system populated on the site 1) obtained from ITM and 
master equations of Redfield form are plotted in Fig.l. 

In Fig.l, we plot the time evolutions of the population 
of site 1 by using the master equations of Redfield form 
and the ITM algorithm as A = A/50, A/5, A, and 5A. 
The initial state is |-0o) = where |1) is the basis state 
of the system. It is shown that the evolutions of the pop- 
ulation of site 1 obtained from the master equations of 
Redfield form are in good agreement with the ones ob- 
tained from the ITM algorithm of 5fc max = 0. We shall 
use <5fc max = 3, and St = 10 fs to investigate the dynam- 
ics under non-Markovian approximation in this paper, so 
in this case, the memory effects of the bath within 30 fs 
are included. In fact, the memory time of the bath with 
Drude-Lorentz density function is longer than this value, 
the memory effects beyond 30 fs will not be included in 
this paper. The time evolutions of the population of site 
1 from ITM of Sk max — 3, are also plotted in Fig.l which 
is quite distinct from the Markovian results. It shows 
that when we investigate this kind of system the non- 
Markovian effects are non-ignorable. In the following, we 
shall use the ITM algorithm with 6k max — 3 to study the 
decoherence and the transfer rate of the EET. The mem- 
ory effects of the bath have not been completely included 
in our calculations but the dominating affects have been 
included, because when we increase £fc max from 3 to 4 the 
results have little change and this change even cannot be 
distinguished by eyes. 



III. TIME EVOLUTIONS OF COHERENCE 
AND EET 

It is reported d, 0| that the decoherence time of the co- 
herent superposition state produced by laser pulses has 
the same time scale of the EET between two adjacent 
chromophores. People thus suppose that coherence may 
help the EET in photosynthesis. On the other hand, the 
theoretical analysis shows that the environment is not a 
hindrance but a help to the EET [13, El- So, we want to 
know further that how the coherence helps the EET and 
how about non-Markovian environment compares to the 
Markovian one in the helps of the EET. In this section 
we shall discuss these kinds of problems through investi- 
gating above dynamical model to mimic photosynthesis 
systems. 

Similar to Sec. II we can easily plot the evolution of the 
coherence, namely, the the evolution of the off-diagonal 
term in the reduced density matrix. It is also convenient 
to plot the evolution of the changing rate of population of 
site 2. The changing rate of the population of the site 2 in 
fact reflects the transfer rate of the EET from the site 1 
to the site 2. Comparing the evolutions of the coherence 
and the transfer rate (denoted by k in this paper) of 
the EET, we can judge whether the transfer rate of the 
EET increase or not as the coherence increase, and as 
the system in a non-Markovian environment instead of 
Markovian one. 

In Figs. 2 and 3 we plot the evolutions of the pi2(t) and 
k with time in different A. Here we set A vary from A/ 5, 
to A/4, A/3, A/2, and e = 100 cm" 1 , A = 100 cm" 1 . 
In this paper we always set 7 = 53.08 enr 1 (7 _1 = 100 
fs), and T = 300 k. The initial state is set in the max- 
imal coherent superposition state |V>o) = + |2)) /\/2 
in Fig. 2, and in the basis state 1-00) = |1) hi Fig. 3. Figs. 4 
and 5 plot the evolutions of the p\2 it) and k with time t 
in different A. Here, we set e = 100 cm -1 , A = 20 cm -1 , 
and A vary from e, to e/2,e/3, and e/4, the initial state 
is |Vo> = (|1> + |2» /y/2 in Fig.4, and |Vo) = |1> in Fig.5. 
In the Figs. 2-5, the time evolutions of the coherence are 
plotted in the upper panels (a) and (c) and transfer rate 
k of the EET are plotted in the lower panels (b) and (d). 
In the left panels (a) and (b) are the results of Marko- 
vian approximation 5fc max = 0, and in right panels (c) 
and (d) are the results of non-Markovian approximation 

max 3. 

The Table 1 is made by using the data of Figs. 2 and 
3, and the Table 2 is made from data of Figs. 4 and 5. In 
the Tables 1 and 2 we numerically show the decoherence 
times and the net area of curve k in different conditions. 
The positive value of transfer rate of the EET means that 
the energy transfers from the site 1 to the site 2, and the 
negative value of the k means that it transfers back to 
the site 1 from the site 2. The net transfer quantity of 
energy from the site 1 to the site 2 is proportional to the 
net area under the curve k (where the net area denote al- 
gebraic sum of the area). From Figs. 2- 5 and Tables 1 and 
2 we can obtain the following results. (1) The transfer 



5k 



4 




t(fs) t(fs) t(fs) t(fs) 



FIG. 2: (Color online) Time evolutions of the coherence [up- 
per panels (a) and (c)] and transfer rate of the EET [lower 
panels (b) and (d)], in Markovian approximation <5fc max = 
[left panels (a) and (b)], and in non-Markovian approximation 
5fc max = 3 [right panels (c) and (d)], for different values of A. 
The initial state is \ipo) = (|1) + |2}) /\/2. The values of other 
parameters are the same as in Fig.l. 



FIG. 4: (Color online) Time evolutions of the coherence [up- 
per panels (a) and (c)] and transfer rate of the EET [lower 
panels (b) and (d)], in Markovian approximation <5fc max = 
[left panels (a) and (b)], and in non-Markovian approximation 
5fc max = 3 [right panels (c) and (d)], for different values of A. 
The initial state is \tpo) = + |2)) /V2. The values of other 
parameters are the same as in Fig.l. 
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FIG. 3: (Color online) Time evolutions of the coherence and 
transfer rate of the EET in Markovian and non-Markovian 
approximations. The parameters in Fig. 3 are similar to them 
in Fig. 2 except the initial state \ipo) = 

rate of energy k from the site 1 to the site 2 are oscilla- 
tions, the coherent term pi2(t) are also oscillations. The 
two kinds of evolution curves almost have the same decay 
times. It means that the system is transferring the en- 
ergy in all of the coherence time. However, the maximum 
of k is nether at the maximal points nor at the minimal 
points of the curve p\2 (t) . This means that the increase 
of coherence does not result in the increase of the trans- 
fer rate of the EET. Thus, we obtain the first conclusion 
in this paper that coherence helps the energy transfer 
through prolonging the transfer time rather than increase 
the transfer rate. As seen in the figures and tables that 
although the time of energy transfer is extended due to 
the increase of coherence its efficiency may decrease be- 
cause the energy will be transferred between the site 1 
and the site 2 back and forth. The breath-like transfer 



FIG. 5: (Color online) Time evolutions of the coherence and 
transfer rate of the EET in Markovian and non-Markovian 
approximations. The parameters in Fig. 5 are similar to them 
in Fig. 4 except the initial state \ipo) = 



of energy in fact is ineffective. (2) From the Tables 1 
and 2 we see that the environment really helps the en- 
ergy transfer. It is clearly shown that an environment 
with stronger damping may stop the breath-like trans- 
fer of energy, which increase the transfer quality of the 
EET. In particular, as the non-Markovian effects are in- 
cluded the net area under the curve k, namely, the trans- 
fer quality of the EET is increased greatly in all cases. It 
means that the environment with non-Markovian effects 
will help the EET. This is the second conclusion in this 
paper. It is worthwhile to point out that as the initial 
state is the coherent superposition state the energy may 
be transferred from site 2 to 1 then the net area of the k 
curve is negative as shown at 410 fs in Table 2. 
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TABLE I: (Obtained from Figs. 2 and 3) Transfer quantity of 
energy and decoherence times in different A as A = e — 100 
cm" . 
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TABLE II: (Obtained from Figs. 4 and 5) Transfer quantity 
of energy and decoherence times in different A as e = 100 
cm -1 , A = e/5. 
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itatively identical to the results obtained from reduced 
hierarchy equation approach in Ref.[12j. The path inte- 
gral method includes the non-Markovian effects of bath. 
By use of the method we furthermore investigated the 
evolutions of the coherence and the transfer rate of the 
EET between two sites. We clearly show that there exist 
the energy transfer in all decoherence time. The increase 
of coherence does not directly result in the increase of the 
transfer rate of the EET. The helps of coherence to the 
energy transfer results from the increase of transfer time 
rather than the transfer rate. The energy transfer be- 
tween two sites within first several hundred femtoseconds 
is like a breathing motion in Markovian environment. We 
have also obtained that the non-Markovian effects play 
important role to the dynamics of the energy transfer in 
the model systems in the initial several hundreds fem- 
toseconds and it results in high efficiency of the EET. 
Our findings from spin-boson model study have impor- 
tant implications for understanding the energy transfer 
in photosynthesis systems. 
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